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Abstract 

We study the relativistic Lee model on static Riemannian manifolds. The model is constructed non- 
perturbatively through its resolvent, which is based on the so-called principal operator and the heat kernel 
techniques. It is shown that making the principal operator well-defined dictates how to renormalize the 
parameters of the model. The renormalization of the parameters are the same in the light front coordinates 
as in the instant form. Moreover, the renormalization of the model on Riemannian manifolds agrees with 
the flat case. The asymptotic behavior of the renormalized principal operator in the large number of bosons 
limit implies that the ground state energy is positive. In 2 + 1 dimensions, the model requires only a mass 
renormalization. We obtain rigorous bounds on the ground state energy for the n-particle sector of 2 + I 
^ . dimensional model. 

B ! 1 Introduction 

Q 

I/*-) ■ The Lee model is a simple field theory model, which requires a mass, coupling constant and wave function 
renormalization [1] . What is so special about the model is that the renormalizations can be carried out nonper- 
turbatively. This is, therefore, a good testing ground for various new ideas and methods on interacting quantum 
field theories. In the original Lee model there are two fermion fields called N and V, assumed to be so heavy 
that their energies are independent of the momentum, and a single relativistic real bosonic field named usually 
as 6. The Lee model is amenable to exact analysis because there are two rather restricting conserved quantities. 
. One of which is the conservation of the total number of the fermion species. Furthermore, the sum of bosons 
and N type fermions is conserved. These highly constrain the theory allowing only a finite number of particles 
£NJ ■ interacting at any given time. If we work with a complex scalar field, the situation changes drastically and the 
model becomes rather difficult [3J. Although the renormalization is performed exactly, it is done so in a small 
number of particle sectors, and it is believed that the same prescriptions will continue to cure the divergences in 
all sectors. This is, of course, plausible since there are no other parameters in the theory. Once the physical V 
£f~) ■ particle is determined as a composite, the coupling constant is determined in such a way to make the scattering 
of N and 9 finite, thereby all the other physical processes should be well-defined (see a source theory approach 
to the model [3])- The model is asymptotically free for d < 4, and this point has been analyzed from the modern 
point of view in [U [5]. There is one subtle point, beyond a certain value of the renormalized coupling there 
appears a ghost state. This problem has been analyzed recently in the virtue of PT symmetry by Bender et al 
in [6] . They show that by an appropriate redefinition of the norm, the ghost state can be turned into a physical 
state. Moreover in [7] , an equivalent hermitian Hamiltonian through a similarity transformation is constructed 
in the context of quasi-hermitian quantum mechanics. Even though the model is defined nonperturbatively, to 
understand the resulting spectrum remains largely as a challenge [8]. 

The model is sufficiently rich, by restricting the total number of fermions to one, we can still get most of 
the interesting features. Moreover, one can assume that these fermions carry no momentum so that they have 
no recoil, and hence assumed to be fixed at the origin. This becomes equivalent to a two states system sitting 
at the origin interacting with relativistic real bosons [HI EH] • This is the version we will be working with so as 
to extend the model to Riemannian manifolds. The nonrelativistic version of this model, which is worked out 
beautifully in the book by Thirring and Henley pj] , is still an interesting case study, yet in this case the coupling 
constant and wave function renormalizations are not needed. The study of the spectrum is still a nontrivial 
problem. To address these issues, there are some attempts in the literature [T5J[T3]. In [H], while looking at 
some nonrelativistic problems which require nontrivial renormalizations, Rajeev introduced a new perspective. 
In this approach one attempts to renormalize the theory by working out the full resolvent in the Fock space 
of the system. The resolvent contains essentially all the information about the model. More interestingly, the 
bound states can be found through the zero eigenvalues of an operator, the so-called principal operator, which 
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is parametrized by the energy in a nonlinear manner. Although one can write the resolvent, it is not possible 
to write down the quantum Hamiltonian of the renormalized theory. In the restricted Lee model, since the 
interaction is at a point, the renormalized model can be considered as a singular extension of the free bosonic 
Hamiltonian. This is analogous to the attractive delta function potential in two dimensions, which requires a 
coupling constant renormalization [15] . There one could also write down the resolvent but not the corresponding 
Hamiltonian. The interaction appears as a kind of boundary condition, this point of view originates from ideas 
of M. G. Krein on operators (see Albeverio and Kutasov [16] for a modern exposition). In [14] . this point of 
view is extended to the nonrelativistic Lee model, inspired from this work we develop the relativistic Lee model 
along the same lines. Having found the principal operator, and thus the resolvent, we can in principle work out 
all the physically important questions for all particle sectors. 

Following the heat kernel based methods, developed in [17], we extend these ideas to the case of manifolds. 
The renormalizations are the same, of course, the resolvent contains information about the geometry through 
the heat kernel. The spectrum of the model is an interesting problem, we only attempt to partially understand 
it for large number of particles and show that the ground state energy remains positive in this limit. Developing 
new approximation methods for estimating the energy levels and scattering amplitudes remains as a challenge. 

The organization of the paper is as follows: In Section ([2]), we will, first, construct the model in flat space- 
time through the approach, introduced in [14j without reviewing it. We show that the principal operator has a 
well-defined limit when the cut-off is removed, and the renormalized operator can be given by the renormalized 
mass and the renormalized coupling constant. Moreover, this limit determines the wave function renormalization 
constant. Afterwards, we specify how to impose the renormalization condition in this approach such that we 
convert the renormalized mass difference into the physical one by fixing the finite arbitrariness, which is left 
after renormalizing the parameters. 

In Section ©, we will apply the ideas, presented in [17] . to the relativistic Lee model on a general static 
Riemannian manifolds. It is shown that the regularization of the ultra-violet divergence in the theory can 
be established through the short-time expansion of the heat kernel if the point interaction is introduced by a 
convolution of the bosonic field with a heat kernel. It is found that the divergence structure of the model in the 
manifold case is exactly the same as the flat case. 

In Section (£[]), we study the asymptotic behaviors of the renormalized principal operator in the large number 
of bosons limit in both flat and manifold cases. In this limit, it is shown that the leading behavior of the theory 
changes substantially. The ultra-static space-time K x H 3 is given as an example. 

In Section (0), we study the model in 2 + 1 dimensions. The advantage is that it is simple and requires only 
a mass renormalization. This allows us to find rigorous bounds on the ground state energy for n-particle sector, 
thus illustrating the power of this method. 

In Appendix, the same techniques are tested in an oblique light-front coordinates. 



2 Relativistic Lee Model in 
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The model which we will construct in this section describes the interaction between a field of relativistic bosons 
and a heavy fermionic source with an internal degree of freedom which actually corresponds to two distinct 
states of the source. Since the source is heavy, we can effectively consider it as sitting at some fixed point in 
space-time. This results in the neglect of recoil for the source, which means that the energies of the states do 
not depend on their momentum [18] . The cut-off Hamiltonian of the model in matrix form is given by 



H t = H \ X + ® X+ + Z(e) X - <8> X-] + #i,e , 



(1) 



where Hq and Hi e are the free and the interaction parts of the cut-off Hamiltonian respectively and are given 

by 

d 3 p 



Hn = 



■^(p)o t (pMp) , 

1 — (7.3 



(2) 
(3) 



At this moment e is an unspecified cut-off prescription, the meaning of which will become clear when we will 
renormalize the parameters of the theory. Here, (f>^ (0) are the positive and negative frequency parts of the 
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real bosonic field, respectively, defined through this cut-off prescription. A more precise dcfiniton of these field 
operators will be given in the manifold case. x± m Eq. ([I]) are the standard spin states, which describe the two 
states of the system. Due to the fact that there can be various divergences hidden inside the theory, we allow 
two states of the system to have different normalizations. With hindsight we choose this to be \- state. 

The theory has a conserved charge, which can be written as, 



Q = -^+/ 7 ^« t (pWp). (4) 



1 - «r 3 , f d 3 p 
(2tt) s 

This means that the theory decouples into independent sectors as ® x+ © •^g™ -1 '' <8> X-- 

The construction of the model is merely based on finding the resolvent of the cut-off Hamiltonian, which 
describes the system completely. While computing the resolvent, one introduces the principal operator &(E), 
that can be regarded as an effective Hamiltonian of the theory. The reason that the Krein formula can be 
applied lies in the observation that there is a constraining conserved quantity, namely, Q. The idea of using 
this operator comes from the fact that the zero eigenvalues of it determine implicitly the bound state energies 
of the theory. The ultra-violet divergence takes place in the theory when the size of the source goes to zero, 
which causes the difference of the energy levels to become infinite. In the following, by renormalization it is 
solely meant to search for a well-defined limit of the cut-off Hamiltonian in matrix form as e — > + . This is 
accomplished by curing the principal operator in the same limit. 

Following Rajeev |14j the cut-off Hamiltonian minus energy is given in a 2 x 2 block form, 



H r — E = 



H -E v ^(i)A(e)4" ) (0) 
yfZ(ej\(e)4 +) (0) Z(e) [H - E + M (e)] 



(5) 



The resolvent is simply the formal inverse of Eq. © and this inverse can be calculated algebraically. If the 
Hamiltonian is parametrized as, 



;;.(Ei , 1.71 



and if the resolvent is parametrized as, 

r m\ — I 

^ 5 ) ' 

then one ends up with the following algebraic equalities, which allow one to calculate the resolvent, 

, -l 



a = a' 1 + a- 1 ^ (d - ba-^Y ha' 1 , (8) 

= - (d - 6a _1 6 t )~ 1 ba- 1 , (9) 

5=(d-ba- 1 tfy 1 =S\ (10) 

$ = d - ba-W . (11) 



Equation (jTTJ) is just the cut-off principal operator and is given by 

..w-^){*-^*)-*"w/^/^' ° (P) 1 ° ,<q) 



r) 3 VMP) ( H o - E) ^Mq) 



(12) 



where E is considered as a complex parameter and the formulae below should be analytically continued to 
their largest domains of analyticity. As one can easily notice the annihilation and the creation operators are 
in the wrong order with respect to normal ordering prescription, so we should normal-order them. After being 
normal-ordered, the principal operator becomes 



1 



2w(p) H„-E + uj(p) 



\-( / f \ °(P) I ri QN 

' " ' (2tt)3 J (2tt) 3 ^/^^Ho-E + w^+u^^/^Jip)} ' 
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The two fractions in the fourth term above can be united by a Feynman parametrization, and then can be 
exponentiated as 



1 



1 



1 



w(p) H - E + w(p) J [(H -E)£ + uj(p)] 



ii Jo 

1 



o n — ±l 



1 - e 



- S (H -E) 



(14) 



In order to evaluate the momentum integral one more identity is needed, and this is the so-called subordination 
identity: 



-sw(p) _ 



du 



1 



-s 2 /4u -uu 2 {p) 



(15) 



Jo " 3/2 

With the help of that identity, we can convert oj(p) in the exponential into ou 2 (p) such that after calculating 
the momentum integral the second term can be given by 



d 3 p 1 



1 



1 



(2tt) 3 2w(p) H -E + uj{p) 4^ 

1 



32tt 2 



u 3/2 J ( 27r )3 



— ttt- / ds se~ s / 4 



^ e^ 2 /°°d SSe - 2 /4. 1 



,3/2 



H - E 

1 _ g-^(Ho-E) 



x _ e -a{H -E) 



V^(H - E) 



(16) 



The use of that identity is to convert the divergence of the momentum integral into a divergence emerging 
from the lower limit of the u-integral, which is now mollified by e, explicitly. The momentum integral is no 
longer divergent, and thus, can safely be computed. If the same calculations are done step by step for the three 
fractions in the fifth term in Eq. (|13p without calculating the momentum integral, one obtains 



OO pCG 



- 7 J=- \ —^ = 11 da I dpi dse -<oM(s+^) e -"(p)W) e -s(H -E) 

^/2U{^ Ho - E + uj(q) + lu{ P ) n Jo Jo Jo 



OO f oo 



— I ds I da (s + a 2 ) I d(3 (s + f3 2 ) I du x 
2tt Jo Jo Jo 



3/2 



du? 



3/2 



(17) 



After plugging Eqs. (jTHJ) and (fTT)) into Eq. (fl"3")l , the cut-off principal operator becomes 

f X 2 (e) f°° e^ 1 "™ 2 f x 2, [l-e-^t^o-B) 

$ e (E) = Z(e) {Ho-E + (M(e) - I du —^r- I ds se~ s ' A 



32tt 2 



,3/2 



y/H(H - E) 



A 2 (e) 
2tt 2 J 



oo poo 



ds / da (s + a 2 ) / d/3 (s + (3 2 ) / du 



3/2 



duo 



3/2 



(18) 



It is now time to renormalize the cut-off principal operator. If the exponential in the second integral term 
is expanded in power series in s, one can notice that the only terms which produce divergence are just the ones 
which are up to order s 2 . On the basis of this expansion, one can redefine the coupling constant and the mass, 
whereby the cut-off principal operator can, easily, be regularized. Thus, we are able to achieve the renormalized 
counterparts of both those parameters and the principal operator. In order to accomplish this, it is appropriate 
to divide the principal operator by the square of the coupling constant. The main difference of the relativistic 
Lee model from the nonrelativistic one resides in not only that there is a coupling constant renormalization 
besides the mass renormalization but also there is a wave function renormalization. Therefore, the ratio of the 
principal operator to the square of the coupling constant, Q(E)/\ 2 , should be renormalized instead of just the 
principal operator, &(E), in the relativistic Lee model. In the light of the above discussion, one can renormalize 
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whole the parameters of the model. The following choices regularize the principal operator by canceling the 
divergences, 



/i(e 



-um poo 



A 2 



and then the principal operator is given by 



. e ^ i~ 2 _ s 2 /4 
A 2 (e) ~ A| ' 32^ J e du u s/2 J Q dss e 

1 1 i /'OO —urn 2 /'OO 

■(e) A|j 64tt 2 7 e u 7 



(19) 
(20) 



<S>e(E) 
A 2 (e) 



Z(e)< 



(JT -£) , Mi? 1 



\ 2 R 32tt 2 



dsse' 3 ' 2 ' 4 



\M#o - E) 



1 - sV^(H -E) + ^s 2 u(H - Ef - e -^(Ho-E) 

^ f- OC /-OO /* OC /■ OO 

— / ds I da (s + a 2 ) / dj3 (s + (3 2 ) I dm 
' 7r Jo Jo Jo Jo 



e -(s+a 2 ) 2 /4ui yoo e -(s+/3 2 ) 2 /4«2 

du 2 - 



3/2 



Jo 



3/2 



d 3 p f d 3 q 



(2tt) 3 7 (2tt) 3 



(21) 



We notice that the subtractions in the second line resemble the regularization of the infinite Frcdholm determi- 
nants. This is analogous to the quantum effective action calculations via regularized determinants in the path 
integral formalism. It is obvious that this operator has a well-defined limit as e — > + when both sides are 
divided by Z(e): 

* £ (£?) _ $ R (E) 



}™ + \ 2 {e)Z(e) " \ 2 R ' 



(22) 



and the renormalized principal operator in terms of the renormalized mass and the renormalized coupling 
constant can be given by 



<S> R (E) (H - E) ^ 1 



A| 32tt 2 



roc ^—um roc ^ 

/ du — — — / dsse~ s ^ 4 - 
Jo u 3 ' 2 Jo 



u 3 ' 2 Jo V^(Ho-E) 
1 - sy/u(H - E) + \s 2 u{Ho - E) 2 - e -«v^(ffo-B) 



— / ds da (s + a 2 ) / d(3 (s + (3 2 ) / du 
2 71 " Jo Jo Jo Jo 



3/2 



^(s+^f/im roo e -(s+0 2 f/Au 2 

i — J du 2 - 



3/2 



(23) 



It is easily seen that this limit also fixes the wave function renormalization constant Z(e) to be equal to A^/A 2 (e). 
Now the elements of the resolvent can also be given in terms of the renormalized principal operator. To see 
this, we look at the expression for a, 



a = -J— + -J— ^->(0) Z(£)A ( V>(0)^— 
Ho-E Ho~E 9t () <f> e (E) 9t y) H -E 



(24) 



If we now use the wave function renormalization constant in this expression, we may take the limit e — > + , 
giving us 



1 



a = 



Ho — E Hn — E 



(0)_%^(+)(0)— - 



Similarly, for the others, we find 



= - 
5 = 



$r(E) 
1 



(0) 



Ho-E 



1 



Hq — E 



(25) 

(26) 
(27) 
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These equations tell us that zero eigenvalues of the renormalized principal operator determine the bound states 
and the corresponding energies as nonlinear eigenvalue equations. Notice that the renormalized operator <&r(E) 
converts a divergent linear problem in the Schrodinger picture into a highly nonlinear but a well-defined problem. 

It is also important to know how the divergences are controlled by the cut-off parameter e in the redefinition 
of the mass and the coupling constant. In order to find it out, the integrals in Eqs. (fT§|) and ([2ll|) should be 
calculated. The cut-off dependent parameters can be given as an asymptotic series in e, 



A 2 (e) 



HR 



1 1 



8tt 3 / 2 



as 



1 



A 2 (e) 



1 

A| 



1 

8^ 



ln< 



as 



0^ 



(28) 



To make contact with the usual perturbative renormalization, we will recast the Hamiltonian into a renor- 
malized part and a counterterm Hamiltonian. We will see that there are no other counterterms needed other 
than the ones existing already in the original Hamiltonian. So as to establish that, we should go back in the 
calculations and replace Z(e) by X R /X 2 (e). The mass term of the source becomes 



Z(eMe) 



1 - cr 3 



-An — 



0-3 



A 2 (e) 



[VR 



(29) 



in which the term Afj, is nothing but the divergent part in the redefinition of the mass. The same replacement 
should also be done in the interaction terms and one can get 



V^)A(e) a + 4-H0) +<t_^ + >(0) 



= A; 



+ 



( - ) (O) + CT --0 (+) (O) 



(30) 



The next step to perform is plug those into the Hamiltonian such that the renormalized Hamiltonian can be 
determined. After plugging them, the Hamiltonian becomes, 



H e — H 
= H 



X+®xi + Z{e) X -®xl +Z(e) A1 (e)^-^ + v ^)A(e) o+^->(0) +<7_#>(0) 



X+ <8> xl + Z(e) X - ® xl + (l*R + X% A/x) + A fl o+^-> (0) + (0) 



(31) 



We know from the theory of renormalization that if one would like to give the Hamiltonian of the theory in 
terms of renormalized parameters instead of bare or cut-off parameters, then the bare Hamiltonian is given by 
the renormalized Hamiltonian containing only the renormalized parameters plus the appropriate counterterms. 
Therefore if we choose the cut-off Hamiltonian as 



H e — Hr + H [Z(e) - l]x- ® X- + AflA M 



1 - cr 3 



(32) 



then the renormalized Hamiltonian of the theory can be given by 



Hr — Hn 



X+&X++X-® X- 



l-Q-3 , . 

— ^ 1- *R 



( - } (O) + ff-0 (+) (O) 



(33) 



The renormalized Hamiltonian Hr should not be confused with what we call the quantum Hamiltonian Hq, 
which determines the time evolution of the quantum system. The resolvent that we have found in the Fock 
space should correspond to the resolvent of the Hamiltonian Hq defined in this Fock space. The existence of 
this Hamiltonian cannot be proved by a straightforward application of the resolvent convergence as is done for 
a different model in [19j . This question is delicate in our case. However its resolvent can be explicitly derived, 
this Hamiltonian may not be written as an explicit formula. 

Although the renormalized parameters had been found, we did not complete the renormalization. Since 
after regularizing the parameters by removing the divergences, there remains a finite arbitrariness :20J. In 
order to fix these finite parts, which results in determining the physical parameters of the theory, one should 
impose the renormalization conditions. In perturbative field theories, these conditions should be imposed on the 
superficially divergent Green's functions to determine the coefficients of the counterterms and one demands that 
Green's functions satisfy them order by order if these conditions are satisfied to lowest order. In our formulation 
we should also specify similar conditions. In this approach the Schrodinger equation is replaced by the equation 
^(E)^ = 0. So a natural choice is the one related to the simple composite which consists of a single boson and 



G 



X+ state giving us a dressed \- state. We can fix the mass difference of \- an d x+- Therefore, we impose the 
following, 

$ R (E = fi p )\0)=0, (34) 

where fi p is the physical mass difference. If the calculations of the principal operator are followed backwards, 
one can obtain a much more compact version of the principal operator. After a little algebra, we get 



$n(E)=H -E + n R 



\2 



d 3 p 



1 



1 



1 | H -E 

w(p) H Q - E + w(p) uj 2 (p) uj 3 (p) 



(35) 



where the dots stand for the normal-ordered interaction term. Incidentally, if we could expand the first term in 
the integral into a power series in Hq — E, the second and the third terms are canceled, leading to a series of 
finite terms. If we add and subtract the second line above with E = p, p and Ho = 0, the resultant operator will 
satisfy the desired condition and fix the finite part of the renormalization. 



\2 
A 7V 



<Z> R (E) =H -E + n p + ^- 



d 3 p 



1 



1 



1 



d 3 p 
(2tt) 3 



-Mp 

w(p) -ftp + w(p) uj 2 (p) w 3 (p) 

1 1 Hq — E] 



$ R (E) = {H -E + fi p )\l + 
d 3 p 



w(p) H - E + u(p) w 2 (p) 
d 3 p 



A 



\2 
A R 

2 
d 3 q 



(2tt) 3 
at(q) 



^ 3 (p) 
1 



1 



w(p)[ffo-S + cj(p)][- Aip +w(p)] w 3 (p) 
1 a(p) 



(2tt) 3 y (2^) 3 ^2^0 H -E + w(q) + w(p) ^2^) 



(36) 



In Section ([U we will use Eq. to analyze the asymptotic limit of the theory with the assistance of the 
asymptotic limit of the principal operator. 



3 Relativist ic Lee Model on Riemannian Manifolds 



First, we will summarize the necessary tools before going into the details of the construction of the model on 
Riemannian manifolds [21] . We consider Ad Riemannian manifold equipped with a metric structure which is 
static. That is to say, there is a timelike Killing vector field and there is a family of spacelike hypersurfaces 
orthogonal to the Killing vector everywhere. Alternatively, there is a coordinate system in which not only are 
the metric components fj fny independent of the time coordinate, but also g 0j = for j =/= 0. 

It is assumed that the action of a bosonic field can be given by 



S = 



d A x 



(37) 



where £ is a dimcnsionless constant and R is the curvature scalar of the manifold. The indefinite analogue of 
the Laplace-Beltrami operator, the so-called wave operator obtained through the covariant derivative, and the 
resulting field equations are given by, respectively, 



□0 = ^ ^y/gdv* 

V9 

U<t> + (to 2 + £R)0 = . 



(38) 
(39) 



Since the metric is static, the field equations can be solved by separation of variables, which results in the 
eigenvalue equation for the operator L, 



L<j) = 500 

L<f>j = uj 2 (j)(f>j , 
where the bosonic field decomposed as 



i=0i(vW*^) + (™ 2 + W 



(40) 
(41) 

(42) 
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The operator L is formally self-adjoint with respect to Cp inner product defined through 



<h) = / <rxy/\g\g m fi(x)M*) 



(43) 



Any function in the Hilbert space defined by that inner product can be expanded in terms of the solutions of 
Eq. (|4"0|) , namely the eigenfunctions of that operator, as 



d/J>(j) (j}(j)(j)j(x) , 



(44) 



where J dfi(j) is the measure and it contains point spectrum or discrete spectrum or both. By means of that 
expansion, the scalar product defined by Eq. (143j) can be given by 



i) 



We have also the orthonormality and the completeness relations. They are, 



(45) 

(46) 
(47) 



The general solution of the field equation can be decomposed into positive and negative parts and they are given 
by, respectively, 



5(j,k) = j ^VMs u >*(x)<fe(z), 



<p(t,x) 
</><+> (a:) 



a(J)4>j (x)e-^ t + a tO>*(x)e^«* 



= (j)j(x)a(j), 



VMi) 



(48) 
(49) 
(50) 



where a(j) and at(j) are the annihilation and the creation operators. A conjugate momentum and a Hamiltonian 
should be defined in order to quantize the field canonically. The conjugate momentum is 



ir(t, x) 



g\do (, 



(51) 



and the Hamiltonian is just the Legendre transform of the Lagrangian. With the help of them, one can calculate 
the equal-time canonical commutation relations both between the field and the conjugate momentum, and then 
between the creation and annihilation operators. 

{cl ) (t,x),ir(t,x')}=i^\g-\g 00 S^(x,xi), [a(j), a?(k)] = S(j, k) . (52) 

The free Hamiltonian in terms of creation and annihilation operators is given by 

Ho = J gWj'MjVOXj)- (53) 

Since the source is heavy and essentially sits at a point in space, one has to find a way to describe this situation. 
We use the same trick which was used in [17] . The interaction is introduced by a convolution of the bosonic 
field with a heat kernel whose index is just a short-time cut-off. In the limit as the cut-off goes to zero, the heat 
kernel becomes a Dirac delta function and hence the convolution in this limit allows us to find the interaction 
occurring at some fixed point in space. Utilizing the short-time behavior of the heat kernel is a nice way to 
analyze and control the high energy behavior of the expressions, so this allows us to deal with the ultra-violet 
divergence in the theory. The cut-off Hamiltonian of the theory on Riemannian manifold, specified previously, 
is 



X+ ®xl + Z{e) X -®xl +Z(e) M (e)i-^ + v /Z(^A( e ) a+4~\w) + <t_^+>(3?) 



(54) 



S 



in which the smeared out positive and negative frequency parts of the field are given by 

4 + \x) = J d 3 g xK e/2 (x,x)^(x), 

4~\x) = f d 3 g xK e/2 (x,x)^-\x), 



(55) 
(56) 



where J d g x = J d 3 x \/\g(x)\g 00 (x), x is a fixed point on the manifold, and e/2 is chosen for convenience. 

Before carrying on, we would like to list the important properties of the heat kernel, [22] . which we will use 
throughout this and the next section: 



K u (x,y) = K u (y,x) , 

d 

LK u (x, y) = —K u (x, y) , 
lim K u (x,y) = S g 3 \x,y), 

d 3 g z K Ul (x, z)K U2 (z, y) = K Ul+U2 (x, y) , 

M 

K u (x, y) > for all u , 
The resolvent is again the formal inverse of the operator. 

_ £ _ //,, 



Symmetry property , 
Heat equation , 
Initial condition , 

Reproducing property , 
Positivity . 



(57) 



(58) 



_y/Z(e)\(e)4 +) (3?) Z(e) [H -E + M (e)]_ 
The cut-off principal operator can be calculated algebraically by the resolvent as in the flat case and is given by 

$ £ (E) = Z(e) |ff - E + M (e) - A 2 (e)0(+> (s)^±-^4~) (5F)} 



= Z{e)lH -E + ii(e)-y{e) / d g xd g y K e/2 (x,x)K e/2 (x,y) 



VMj) VMM 



0i(^)</ , fe(2/)«(j) 



Hn-E 



at(fc) 



(59) 



Henceforth, the same game is played in order to renormalize the theory. First of all, one should normal-order 
this object by letting the creation operator stand on the right and the annihilation operator stand on the left 
in the fourth term in Eq. (|5"9")l . If the following operator equalities are used, 



a(j)- 



H -E 
1 



a?(k) = a 1 (k) 



1 



H - E + oj(k) ' 
1 



H -E + u(k) H - E + w(k) + u(j) 
then the principal operator becomes 



a (j) , 



(60) 
(61) 



$ e (E) = Z(e)^H Q -E + v(e)-\ 2 (e) J d 3 g xd 3 g y K e/2 (x,x)K e/2 (x,y) J 



H - E + Lu(k) + u(j) 



d^{j)<t>j{x)4>* j {y) 



1 



1 



' > J 2u(j)Ho-E-u(j) 
(62) 



We will again use a Feynman paramctrization and do an exponentiation. After that we compute the Feynman 
integral as 

f i 



1 



1 



1 



uj(j) H -E + Jo ~* [(H - E)C + u{j)f 

f>l POO 

dC / dsse- s ^e~ s ^- E ^ 
Jo 

I 



dse- sujU) - 



H -E 



I _ e s(H -E) 



(63) 
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By means of the subordination identity, can be turned into uj 2 (j) which allows us to convert e sw ^ into 
a heat kernel via sandwiching it with the eigenfunctions of the operator L: 



K u (y,x) = J d f i(j)<f> j (x)<f>*(y)e- uu «. 
The reproducing identity also allows us to combine the convoluted heat kernels as 

K u+e (x,x) = Jd 3 g xd 3 g yK e/2 (x,x)K u (x,y)K e/2 (y,x). 
If all of them are taken into account, we reach the following form of the principal operator, 



$ e (E) = Z(e)iH -E + »(e) - X 2 (e) J d g xd g yK e/2 (x,x)K e/2 (x,y) 
dfi(j) dfi(k) I 



4>j( x )<f>k(y) a {fy 



A 2 (e) 



OO /-00 



du / dsse s / i K u+e ( 



H -E + u(k) + u){j) 

I _ e -s^l(H -E) 



X. X 



V^(H - E) 



a(j)=VMj)J dlzcj>*{z)^ + \z), 
a\k) = ^/Mk) I ' dlz cj> k {z)^-\z) . 



(64) 



(65) 



(66) 



40? Jo 

We can also exponentiate the fraction in the fourth term: 

1 r°° 

- — T -= / dse~ s ^ k h- s ^e-^ H °- E \ (67) 

H - E + w{k) + uj(j) J 

Moreover both a(j) and a)(j) can be given by the field itself by an inverse transform: 

(68) 
(69) 

Equations (I67p . (I68p and applying reproducing property one more time brings the principal operator to the 
following form 



$ e (E) = Z(e)yi -E + 11(e)- 
X 2 (e) 



^ ^ / du I dsse S ^ 4 K U+ Jx,x) 

Jo Jo 



V^(Ho - E) 



-in 



d 3 xd 3 y / ds s 2 I dui 



f „'/4«i poo e _ s ^/4u 2 
— 



3/2 



X ^e/2+ Ml y)^e/2+n 2 z)</> ( } (?/) 



(70) 



We are now ready to renormalize the cut-off principal operator via redefining the cut-off dependent parame- 
ters in terms of renormalized ones and divergent parts so as to cancel the divergences emerging from the normal 
ordering. As in the flat case we, first, determine which powers of u in the it-integral produce divergence in the 
fourth term. In order to do that one should use the short-time expansion of the heat kernel which generates 
some inverse powers of u. We combine it with the powers of u coming from the expansion of the exponential in 
s. Well-known short-time expansion of the heat kernel (see |22j . for example) is given by 



K u (x,x) 



U nu )d/2 zZ 

v ' n—0 



a n (x)u n , 



(71) 



where a n are universal polynomials in the curvature tensor, its covariant derivatives and various contractions 
thereof. The mass term in laplacian does not affect the asymptotic expansion. It can immediately be seen that 
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only the first term in this short-time expansion contributes to divergences when it is combined with the factors 
coming from the exponential. The following choices are sufficient to kill all the divergences, 







1 , 


A 2 (e) 




J 






1 f 00 




= A| + 


2 7 


1 


1 




A2(e) 


A R 






1 


- / 




A R 





dss 2 e s / 4 K u+e (x, x) 



duK u+e {x, x) , 

/>OC />00 

' duy/u / dss 3 e _s "/ 4 i ; i' tl+e (s,x) 



(72) 



(73) 



By taking the asymptotic limits of the integrals above, one can find out how the cut-off parameter e controls 
the divergences and one gets 



A*(e) 



HR 
A R 



1 1 



8tt 3 / 2 
1 



A*(e) " 



8tt 2 



ln(e) 



a.s 



as 



>0 4 
0+ 



(74) 
(75) 



It is striking that these are exactly the same results which we have found in the flat case. This equality arises 
from the fact that in this language short-time behavior captures the high energy behavior. Therefore, we would 
not expect any contribution from the geometry itself, only extreme local structure which is Euclidean determines 
the divergence. It is also easy to conclude that point from the short-time expansion of the heat kernel since 
only the first term contributes to divergences. Moreover, the first expansion coefficient ao does not contain the 
curvature scalar, it is just equal to 1. Yet, the geometry is very important for the spectrum of the theory. The 
principal operator is given in terms of the heat kernel at arbitrary times as well as its values at separate points. 

After replacing the parameters by their renormalized counterparts and successively taking the limit e — > + , 
we can obtain the renormalized principal operator, which is given by 



$ R (E) _ (H - E) m 1 



A l 



— tt j= / du dsse~ s ^ 4 K u (x,x) 

A R 4 V T Jo JO 



1 



V^(H - E) 



[1 - s^/u~(H -E) + -s 2 u{H - E) 2 - , 



-s^(H -E)] 



ds , 



du\ 



2 /Am 



3/2 



duo 



3/2 



x K Ul (x, y)K U2 (x, x)(j) { ) (y)e~ 



s(H -E] 



4> {+) {x). 



(76) 



After imposing <&r(E = fJ> p )\0) = and doing little algebra, one can also obtain the principal operator in terms 
of physical mass difference as in the flat case, 



$ R (E) = (H -E + fi p )\l 



\2 
A R 



1 



>(j) 3 



>{j) [H -E + Lj(j)] [-Mp + w(j')] 



d^(j) dfi(k) 



^i{x)(t>*k{x)a) (k)— — 7-t — - 

VM7) VMfc) 'H -E + cj(k)+u(j) 



(77) 



We see that the renormalized Hamiltonian, after the same calculations done in the flat case, is given by 



H, 



Ho 



X+®X+ + X- 



fJ-R- 



1-0-3 



+ A/ 



(78) 



4 Asymptotic Limits 

In this section, we will study the asymptotic behavior of the operator $>r(E) in the limit of large number of 
bosons, n — > 00, and the flat case is our starting point. Combining the fractions in Eq. (|36[) by Feynman 
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parametrization, exponentiating the resultant fractions, applying the subordination identity to e SUJ ^ and 
taking the momentum integral, successively, brings the principal operator to the following form 



I A 2 f 1 fi~£ f°° i f°° 
* R (E) = (H -E + ii p )<l + 2^J J d( J ds -J du- 



-l/iu-us 2 



-s[(H -E)t-n P C] _ 1 



(79) 



The u- integral is just the integral representation of the modified Bessel function of the second kind K 2 (ms) 
multiplied with 8m 2 s 2 . After letting u — ► us 2 , calculating the u-integral and scaling sass^ s/m, one gets 



d( / dssK 2 (s) 
Jo 



e -a[(H -E)(-n P C]/m _ 2 



(80) 



We can add and subtract e _s to the right hand side of this expression, which regularizes the s-integral and 
also generates a constant — C, which is numerically computable and approximately equal to —2.67. Thus, the 
principal operator is 



dt 



d( / dssK 2 (s) 
Jo 



e -s[(H -E)t-fi p £]/ m _ e - 



(81) 



Being calculated by some mathematical software, the s-integral is equal to 



bam 



8a: 



+ (-12 + y/2 - 2a, /m) m 2 - V2to 3 / 2 V '-a + to 



20a/1 — a/m(m — a) [m(a + m)] 3 ^ 2 
32(a - ro)(a + m) 2 V -a 2 + m 2 3 F 2 (l, 1, 5; 2, |; ^p) 



20^/l — a/m(m — a) [m(a + to)] 



3/2 



(82) 



where a = (H a — E)£ — fi p £. After converting the Hypergeometric function into an elementary function and 
doing some simplifications, the principal operator can be given by 



^ R{ E) = {H,-E + ^){l-C^ + ^j\^^ 'd( 



TO 2 



(C-C) (C-C-2m) 



+ 



(<-< 



(<-<)( 



2to - C + C 



(c-c) 3/2 (2to-C + c) 



~\ 3/2 



7T- 2arcsinl/l - -4- + -4- 
V 2to 2to 



(H -E)£-m 



(83) 



where to = — and C _ 

Whether £ is between the limits of the ^-integral or not is important for calculating this integral. Taking 
the integration interval of the ^-integral and H > m into account tells us that £ is in the integration interval. 
This could cause poles since the denominators has some powers of ( — (. If this is the case, then the integral 
should be defined either by a principal value prescription or by a Hadamard finite part prescription. In order to 
answer this question, one should expand the integrand around ( — ( in series. For this expansion, logarithmic 
form of the inverse trigonometric function is more suitable: 



2 arcsin \ / 1 



2to 2m 



= 2i In 




(84) 



The scries expansion of the combination of the first and the second term is 

to 2 to 23 C-C 



lim \ 



< 3 (C-C) (C-C-2m) 2 (C - C) 



12 



8to 



+ 



o(c-ey 



(85) 
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and the expansion of the third term is 



lim 



2m - ( + C 



c-c 



~\ 3/2 



2m - C + C 



~\ 3/2 




~\ 3/2 



(86) 



It is astonishing that not only the singular parts but also the constant parts of the integrand in the expansion 
cancel each other and that limit of the integrand is just given by 



lim (integrand) = (( - C) + 0(C - C) 3/2 . 



(87) 



Although £ is between the integration limits, series expansion tells us that the £ integral is an ordinary integral 
since the integrand does not have any poles at £ ~ £. Thus, we do not need to introduce any prescription to 
compute this integral. After tedious calculations, the exact principal operator can be obtained as 



A 2 



7\ A 



4tt 2 



A 2 « 



E - m) (Hq - E + m) In 



H - E-m 



2m 



H - E + m 



2m 



m — ii p 
2m 



A 2 h 



A « f (q) 



«(P) 



d 3 p 

J (2^) 3 ^2ZJ(^)^o-^ + ^(q)+w(p) VMp) 



In the flat case, the asymptotic behavior of the principal operator in the limit of large number of bosons, 
that is Hq > nm S> m > fj, p , is, then, given by 



<Z> R (E) ~ /r 



1 - ( C- 7 - +ln2 



\ 2 

47T 



A 



K H n \n 



47r' 



Ho 



(the normal-ordered interaction term) 
(the lower order terms in Hq) . 



(89) 



This asymptotic behavior has a striking feature, the interaction term is positive, multiplied by a minus sign gives 
a negative contribution, and the leading term of the renormalized principal operator is also negative. Whatever 
the leading behavior of this interaction term is, these two terms are enhancing the negative value of &r(E). We 
can show the positivity of the interaction term in general by studying the same term in the manifold case. The 
interaction term in Eq. (|76p can be written as 



4tt 2 



POO 




/ dss 2 


J d 3 g y J dui 


IQ 








x e - s{H °- 


E) 









2 /4«i 



3/2 



K Ul (x,y)4> {+ \v) 



t 



du 2 



-s 2 /4« 2 



3/2 



K U2 (x,x)^+\x) 



\ 2 
A R 

At: 2 



dss 2 A\s)e- s ( H °- E ^A{s) 



(90) 



>o 



Since the integrand is positive, the interaction term is positive-definite apart from the minus sign in front. This, 
in turn, implies the operator to have a negative-definite sign. Henceforth, the operator &r(E) can not have 
zero eigenvalues for E positive but small. For large number of particles this proves the positivity of the energy, 
which is extremely important for the stability of the theory. 
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Secondly, we will analyze the behavior of the principal operator on a general ultra-static manifold in the 
same limit. Having done similar calculations, Eq. (|77j) becomes ready to be studied in the limit n — > oo. 



<£ R (E) = (H Q -E + Li p )\l-C(x,m)-^= , 

4y7r 4y7r Jo 



Ar I dd dC I dss 3 



(I 







- I du -j- K u (x, x) 

M J/ 



e -a[{H -E)i-nQ _ e -«ra 



(91) 



Appropriate scalings of the variables in the above equation allow us to take that limit and the operator is given 

by 

f \2 \2 pi pl-(, poo 3 

$ R (E) = (H -E + » p )\l-C(x,m)-^ + -!j= / / dC / <*« 7 « 



-um 2 j (nm) 2 — s 2 / '4u 



< I du- 

/Q 



.3/2 



-A" 



u/(r, 



! (x, x) 



e —s[(Ho—E)$—fj, p Q/nm _ ^—sm/nm 



(92) 



The asymptotic behavior of the heat kernel is given by 



lim K, 



uj (nm) : 



: (x, x) 



(nm) 3 



(4ttw) 3 / 2 ' 

Plugging the equation above into Eq. (|92]) allows us to take the u-integral and we get 



{a 2 a 2 r 1 r 1 ^^ r°° i /s\ 



e -s[(H -E)i-iJ, p C]/nrn _ g-s/n 



(93) 



(94) 



We should, now, be careful about the asymptotic expansion of the integral. Although asymptotic behavior of 
the function K2{s/n) can be used for s small enough, we are not allowed to use it when s becomes comparable 
with n because the other multiplying factors do not decay sufficiently fast with s. Since the upper limit of the 
s-integral is at infinity, this is the case. However, if we rescale s with n, this integral takes a form which is 
independent of n. Therefore, this expression becomes the same expression which we have found already in the 
previous case whose constant term C is, basically, replaced by tt 3 / 2 C(x, m). If one takes the next term in the 
short-time expansion of the heat kernel into account, then it can be seen that the contribution coming from that 
term is of the order of 1/n 2 , which is much smaller and hence neglected. Yet, there comes a contribution from 
the expansion of the first exponential, which results in a new constant C", multiplying Hq. Thus, the leading 
behavior of the operator <&#(/?) in the asymptotic limit Hq » ra on a general ultra-static Ricmannian manifold 
can be given by 



1 - 



A 



7 



, ., . .. C(x,m) + C - - + ln2 

47T 2 \ 3 

— (the normal-ordered interaction term) 
+ (the lower order terms in Hq) . 



X R. Hq \ n (El 

47T 2 V 1TI 



(95) 



At this stage, we are unable to give precise asymptotic analysis of the normal-ordered interaction term, which 
requires a delicate study. We would like to call readers' attention to the fact that the same remarks, which 
have been done in the flat case, are also valid for the relativistic Lee model defined on a general ultra-static 
Riemannian manifold. 

At last, the manifold defined as Ai = IxB 3 will be considered as an example. H 3 is, here, a three 
dimensional hyperbolic space. The reason why we study this manifold is based on the fact that its heat kernel 
is one of the simplest and explicitly known heat kernels. 
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The heat kernel of the hyperbolic space H™, found in [23], and the diagonal heat kernel of H 3 takes the form 

1 



K u (x,x) 



l ila P&lMi c -a 2 u-p( S ,v) 2 /4u 



(Airu) 3 / 2 v^x sinhp(x,2/) 

e -a 2 u 



(4tto) 3 / 2 ' 



(96) 



where p(x,y) = dist(x,y) is the geodesic distance on H 3 and —a 2 is the constant sectional curvature. Having 
used the diagonal heat kernel in Eq. (|9"Tj) . the following operator could be obtained, 



* R (E) = (H -E + fx p )il 



\2 

32tt 2 



dss" 



du ■ 



2 /4u 



(97) 



This is the same result, which was found already in the flat case with m 2 replaced by m 2 + a 2 . Thus, the 
space H 3 modifies the mass term of the exact principal operator only. After this slight modification, the exact 
principal operator takes the form 



$ R (E) = (H q -E + tip) 
x In 



C-- 



7\ A 



3 J 4tt 2 



2^ J (h -E- Vm 2 + a 2 ) (h -E + ^Jm 2 + a< 



Hq — E — \Jm 2 + a 2 H - E + \/m 2 + a 2 



2\/m 2 + a 2 



2\/m 2 + a 2 



I \j (V m 2 +a 2 - pp^j (V m 2 + a 2 + arccos \j 



\/m 2 + i 



1\Jrn 2 + a 2 



(the normal-ordered interaction term) . 



(98) 



It is easy to see the asymptotic behavior of the operator $>r(E) a simple modification is sufficient to calculate 
it. Hence, one get the following, 



$ R (E) ~ H 



A 



- -In2-C 



4tt 2 V3 

(the normal-ordered interaction term) 
(the lower order terms in Hq) , 



H 



V to 2 + a? 



(99) 



and the normal-ordered interaction term, of course, changes drastically (see Eq. ijTB"])). 



5 Lee model on 2 + 1 dimensional Riemannian manifolds 



In this section, we make a digression to an analysis of the two dimensional version of the Lee model. Our 
purpose here is two fold, we first would like to show that the two dimensional model is much simpler, which 
only requires a mass renormalization and secondly we would like to illustrate the power of this approach by 
obtaining an explicit bound on the ground state energy in each sector. 

We write the model on a Riemannian manifold in the matrix form by using a heat kernel cut-off function: 



He — E 



H a -E 

,(+) f-TS 



\4~\x) 



XcPt'ix) [H -E + fi(e)] 



(100) 



The model now neither requires a coupling constant renormalization nor a wave function one. We take the 
resolvent in the same way as before and find the principal operator as, 



$ e (E) = i H -E + 11(e) - X 2 / dt g xd*yK e/2 (x,x)K e/2 (x,y) 



-jm 7m Uxmy)aHk) ho-e+ i( k ) + w) 



dfj,(j)<pj(x)<f>j(y) 



1 



1 



2w(j) Ho-E- Lo{j) 

(101) 



15 



Following the same steps in the 3 + 1 dimensional case, we end up with, 



( \2 roo roc 

$ e (E) = I H - E + n(e) - —= du dsse" 8 / 4 K u+e (x,x) 

I 4^7T J Jo 



47T 



d 3 g xd 3 g y / dss 2 du\ ^ — / du 2 

Jo Jo Ui Jo u 



V^(Ho - E) 

2 /4u 2 



3/2 
2 



x K e/2+Ul (x, y)K t/2+U2 (x, x)^ (y)e- s ^- E U (+) (x) 



(102) 



Using the behaviour of the heat kernel on a two dimensional Ricmannian manifold, 

1 

K u {x,x) ~ — — rVa„(i)ti", 
(4™) ^ 

we see that the principal operator becomes finite if we define a mass renormalization given by 



/"(e) = Mi? + 



A 2 



2 t'OG 



I du I ds s 2 e~ s2/i K u+e (x, 
Jo Jo 



a r 00 



As a result we find the renormalizcd principal operator as 

V2 



\2 /■oo /■oo 

$«(£) = (H -E + ii R ) -—^\ du ds se~ s /4 K u (x,x) 
x [1 - sV^(^o — E) — e -^(H -E)} 



r too roc e -s 2 /4«i /■cxi 

-^J< X < y Jo """Jo dUl —Jo 
x K Ul (x,y)K U2 (x,x)^-Hy)e- s ^- E U i+) (x) . 



V^(Ho - E) 



-s /4m 2 



3/2 



(103) 



(104) 



(105) 



If we now impose the physical mass condition $(E = fj, p )\0 >— 0, written in the eigenfunction expansion, we 
end up with, 



$ R (E) = (Ho-E + ,i p ) [l + y / M.l) PiWhtx) 
dfj,(j) du{k) 



- A 



<j>j{x)4>l(x)a) (k) 



Lj(j)[H -E + u(j)} [-n P + w{j)} 
1 



H -E + u{k) + 



(106) 



The change in the renormalized part is important, if we recall that u p < this part is actually always positive 
for E < nm (the interesting case from the bound state spectrum point of view) . Thus the interaction term now 
competes with these two terms. If we evaluate the answer for the flat case we see that it is given by 



Z R (E) = (Ho-E + M + f In \ H o-E + m 

47T m — ii p 



d 2 p f d 2 q at(q) 



a(p) 



(2tt) 2 J (2tt) 2 ^J{q) Ho - E + Lu{q) + u(p) ^J{p) 



(107) 



Since the flat case is sufficiently important we will give a bound on the ground state energy for all particle 
sectors first and discuss the general case of manifolds afterwards. Note that if we can show that the principal 
operator becomes positive for sufficiently small values of E, this means that it is invertible, hence, it cannot 
have a zero eigenvalue beyond that value. This give us a lower bound on the ground state energy. To accomplish 
this we rewrite the principal operator in the form, 



Q R (E) = K(E) - U(E) . 



(108) 
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where 



K(E) = (H -E + Mp ) + — In 

47T 



H - E + m 



Hp 



U(E) = A 2 



d 2 q a^(q) 



a(p) 



d 2 p 

W? J (27r) 2 ^2^#o-£ + ^(q)+cj(p)^2^ 



(109) 



Note that for real values of E, we can drop the logarithm and the resulting operator is smaller than K{E). 
Thus following Rajeev |14j . we write an inequality of the form, 



&r(E) > K(E) - U{E) = K(E) 1/2 (1 - K{EY' rz U{E)K(EY' rz )K{E) 



-1/2, 



-i/2\ 



a/2 



(110) 



where K{E) — Ho+fi p — E. Hence to show that the operator Qr(E) to be invertible, it is sufficient to impose the 
condition, \ \U(E)\ \ = \\K{E)- 1 I 2 U{E)K(E)- 1 I 2 \\ < 1. This will impose a condition on the ground state energy. 
If we write this out explicitly, after commuting the square root operators with the creation and annihilation 
operators of the interaction term, 



U{E) = A 2 



d 2 p f d 2 q 



(2n) 



(2tt) 2 



at(q) 



o(p) 



y^A^j [H -E + fi p + w(q)]i/3[iT - E + w(q) + u(p)][H - E + Mp + w(p)]V2 



(111) 



Now we can use the inequality Hq > (n — l)m in the n boson sector inside the operator and replacement of it 
results in a bigger operator function. Call this \ — (n — l)m + fx p — E and for n > 1 we find as a result, 



U(E) < A 2 



d 2 p f d 2 q 



(2tt) 2 J (2?r) 



ot(q) 



a(p) 



^(qT [X + ^(q)] 1/2 [x - M P + "(q) + W(P)][X + ^(P)] 1/2 v/Mp) ' 
If we now use an extension of the Cauchy-Schwartz inequality to the Fock-Space operators, we find 



\\U(E)\\<^n\ 2 



d 2 p f d 2 q 



(2tt) 2 J (2tt) 2 w(q)[x + w(q)][x " M P + ^(q) + ^(p)] 2 [x + ^(p)Mp) 



1/2 



(112) 



(113) 



We now note that i/p 2 + m 2 > |p| = p and m > /j, p , and replace some the terms by these lower ones and thus 
preserving direction of the inequalities, 



\\U(.E)\\<\n\ 2 



pdpdflp f qdqdfl q 



(2tt) 2 J (2tt) 2 pq[x + q + p] 2 [x + q][x + p] 

Let us scale the momenta by p = xVi 1 — XQ we nnd 



1/2 



C30 /-OO 



dpdq 



1/2 



/o Vo [l + g + p] 2 [l + ?][l+p]. 
And the last integral is finite, let us call its value as C, we then impose the condition, 



(114) 



(115) 



< 1 



8n 2 X 

which guaranties that the < 1 This implies the rigorous inequality on the ground state energy, 

A 2 nC 



E gr (n) > {n — l)m + fj, p — 



8ir 2 



(116) 



(117) 



If we want the energy to be positive in all sectors this in turn brings about a bound on the coupling constant. 
In fact, for a global stability we should have the energy to be bounded by (n — 2)m. But the present analysis 
is too crude to get a bound of this form, that requires a much more delicate analysis. 
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Next we will work out the same problem for the Riemannian manifolds, it is simpler to work on the eigen- 
function expansions. We follow the same approach and estimate the leading behavior of the term resulting from 
the renormalization. 

The denominator of the second term in Eq. (|106p can be united by Feynman parametrization as, 

* iw . Jn „ r« fMf, feg 



2 J V / wCj)[ifo-^ + wO')][-MP+wC?)] 2 Jo Jo J [uj{j)(H - E)£ - fx p C] 

After converting the fraction into an exponential, utilizing subordination identity and the definition of the heat 
kernel, this term becomes, 

\2 pi poo poo -s 2 /4u~-um 2 ~s(H -E)£+ f i p sC 

d£ dC dss 3 / du ^j- K u (x,x). (119) 



4^ Jo Jo Jo Jo u3/2 

Let s — » s/(nm) and u — ► a/(iim) 2 , we obtain, 

d£ / rfC / dss 3 / jj^ K u/{nm) 2(x,x) . (120) 



The asymptotic behavior of the heat kernel for large n is given by, 



lim K u/ , nm y(x,x) ~ — . (121) 



After pluging this asymptotic behavior, one gets, 

a 2 r 1 „ r 1 ^ 



e 



^ / d C / / du — . (122) 



lQTT 3 / 2 nm J Q J Q J Jo u 5 / 2 

Performing the integrals results in, 

~H -E + m~ 



A 2 1 ln 

4tt (H - E + fi P ) 



(123) 



Taking the overall factor (Hq — E + fi p ) into account, we find the same answer as the one in Eq. (|107|) and 
the leading contribution of the renormalization to the principal operator in the large number of particles limit 
results in 



A 2 

Q R (E)~ (H -E + ti P ) + —\n 



Ho - E + m 



m — 



,2 f d Ki) d K k ) j. /_ u «« ui\ 1 r-\ 

-A / — . — . (hAx)<h k (x a 1 (fcj— ttt Trza(j) 

J ^/Mjj VMfc) ; H -s + ^(fc) + ^(j) u; 

+ (smaller order terms) . (124) 

The term in Eq. (|118|) is always positive, and we see that its leading term is of smaller order. Hence can be 
dropped out safely without affecting the inequalities, 

$ R (E) >{H Q -E + m p ) 1/2 [1 - U(E)](H -E + tip) 1 ' 2 . (125) 

We will work on a noncompact manifold, for the compact manifold case the zero mode should be worked 
out separately. We expand the U(E) in the eigenfunction basis, 

U(E) = A 2 J dti(j)d t i(k) 

j2^[H -E + ti p +u{j)]V*[Ho-E + u{fi+u{k)^^^ 1 ' 

Here we have oj(j) = ^/ct 2 + m 2 and we introduce again \ = ( n — l) m ~~ ( we drop the fi p for simplicity). 
Recall that Ho > (n — l)m, we use this inequality, and the Cauchy-Schwartz inequality to find, 

A 2 n 



ll^)ll< 2 



+ w(j)][x + + ^( fc )] 2 [x + w(fe)]w(fe) 



(127) 
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We use the following crude inequality, 

(x + + ^( fc )) 2 > (x + w(i))(x + ^( fc )) » 
which implies the opposite inequality for the inverse, 



(128) 



W(E)\\< 



< 



X 2 n 
X 2 n 



dfj,{j)d(j,(k) 



-,1/2 



u(j)[x + u(j)} 2 [ X + uj(k)} 2 u(k) 



\M*)\ 2 



w(j)[x + ^(j)] 2 



We now employ a Feynmann parametrization, 

1 



2(d( 



^(j)[x+^(j)] 2 y [ W (j)(i-o + (x+^(j))c] 3 



(129) 



(130) 



To make contact with the heat kernel we employ an exponentiation and then use the subordination identity to 
get, 



\\U(E)\\< 



1 



d(( / dsi c 



du 



e -s /4u-ui 0")u 

d ^') 



-l^(3?)l< 



Recognizing the heat kernel as, 



K a (x,x) = I d^ 3 )\m\ 2 ^ U)u 

we can rewrite the desired inequality as, 



1 



\\U(E)\\< — dCC s ds du[e~ m u K u (x,x)} 



-s 2 /4u 
«3/ 2 ' 



We note that for Cartan-Hadamard manifolds there is a nice inequality for the heat kernel [23] . 

K u (x,x) < — , 
u 

where C\ is a positive constant related to the geometry. This in turn implies for these manifolds that 



\mm\\< 



da 



ds / due' 



2 u d e -s 2 /4« 



.3/2 



If we drop e m u term the integral can be easily found, we scale the u variable as sv and find, 

\ 2 n 



\\U(E)\\< 



1 pi f>OC f>OC (-1 



x 



(131) 
(132) 

(133) 
(134) 
(135) 

(136) 



If we impose the condition, Ci^ 11 < 1, then we have no zeros for <&r(E), and this implies a bound on the 
ground state energy, 



E gr (n) > (n ~ l)m — C\\ 2 n . 



(137) 



This shows that there is a rigorous bound on the ground state energy of the n particle system. Again, one 
expects that these bounds are weak, that is, a better physical approximation should prove a better bound. 
Nevertheless the bounds that we found illustrate the power of this approach clearly. 



6 Conclusion 

In this paper, the construction of the relativistic Lee model on static Riemannian manifolds is studied. This 
construction is, basically, based on introducing an operator, the so-called principal operator, and renormalizing 
it successively [14) . Moreover, it allows us to renormalize the theory nonperturbatively. This operator, which 
can be regarded as a kind of effective Hamiltonian of the theory, converts a divergent linear problem in the 
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Schrodinger picture into a highly nonlinear but a well-defined problem. Since it is found through the resolvent 
in the Fock space, it is valid for all particle sectors of theory. Analysis of the behavior of the principal operator 
in different regimes can allow us to obtain definite information about the spectrum of the theory since the zero 
eigenvalues of the renormalized operator implicitly determines the bound state energies. Renormalization in 
this construction is established in two stages. First stage is identifying the divergences in the theory, which are 
tamed by a cut-off at the beginning, and then curing them by redefinitions of the appropriate parameters of the 
model. We show that the principal operator is free of divergences when the cut-off is removed. The second stage 
is specifying the renormalization conditions since there remains a finite arbitrariness in the definitions of the 
renormalized quantities after regularization. Since the renormalized mass of the source /ir should, intuitively, 
be related to the physical mass at the lowest number of particles sector, we believe that a natural choice is to 
impose this condition on the renormalized principal operator. So we choose /i p as the lowest energy solution of 
the equation &r(E)\0) = and replace by this physical parameter. 

As shown, renormalization in the manifold case is much more complicated than the one in the flat case. 
The ultra-violet divergence in the theory is identified through the short-time singularity of the heat kernel, 
the short-time expansion of the heat kernel allows us to determine how to renormalize the bare parameters. 
Only the first term in the short-time expansion contributes the divergences and these can be absorbed in the 
redefinitions of mass and coupling constant. As known, mass and coupling constant renormalizations are not 
sufficient to let the theory be free of divergences so a wave function renormalization is needed. To fix the wave 
function renormalization constant, we start with a Hamiltonian in which a different normalization of two states 
of the system is allowed. In that way, we do not need to change the normalizations of the spin states after 
renormalization. The well-defined limit of a suitable combination of the cut-off dependent principal operator, 
coupling constant and wave function renormalization constant dictates the form of the constant Z(e). The 
divergence structure in the manifold case is the same as the one in the flat case. This is, actually, not a 
surprising result and it stems from the fact the divergence in the theory is an ultra-violet type. We also analyze 
the model in an oblique light-front coordinate system as a case study in Appendix. Same results are obtained, 
which encourages us to confirm the results found in [21]. 

There is another unconventional alternative; where we set the wave function renormalization constant Z(e) 
to — A|,/A 2 (e). This will make Z(e) positive below a certain value of the cut-off e, hence the lower block of the 
Hamiltonian multiplied by a positive divergent number. It will change the off-diagonal blocks into operators 
multiplied by an extra i. To make the Hamiltonian hermitian on C 2 •T'b(W), we should define it through the 
operator, 

H e — E = \ . ¥° ~ E (+) V^JA^O) 1 

It is an intresting alternative to study. 

In Section (jl]), we calculate, first, the exact principal operator in the flat case, and then analyze the asymp- 
totic behavior of it in the large number of bosons limit. The analysis shows that the renormalization process 
changes the leading term distinctively with respect to the free Hamiltonian and it takes the form —HqIuHq. 
This seems to change the dynamics of the model drastically. Therefore one should be very careful how to define 
the quantum Hamiltonian from the constructed resolvent. Another astonishing characteristic of this result is 
the sign of the leading term, which is negative. Since the normal ordered interaction term has also a negative- 
definite sign, the total operator is negative-definite. This implies that the ground state energy is positive. In [25] 
it is shown that the quantum effective action of the large- N Yukawa theory also takes a similar multiplicative 
contribution to the kinetic term. We, therefore, believe these results call our attentions to the point that the 
quantum field theoretical models should be examined in much more detail at the functional level. 

In Section ([5j, to show the power of this approach, we look at the 2 + 1 dimensional model, which only 
requires a mass renormalization and simpler. The model seems to have no ghosts. The cut-off Hamiltonian is 
well-defined. The renormalized resolvet allows us to give a rigorous bound on the ground state. The existance 
of the quantum Hamiltonian can be proved by the methods in Ref. [Hj in 2 + I dimensions. 

How to construct the relativistic Lee model on a general static Riemannian manifold is addressed so far in 
this paper. However, the present analysis does not give adequate information how the spectrum of the theory 
can be build up. Although naive scaling arguments for the normal-ordered interaction term suggest that it 
gives a contribution of order n, a scrutiny of this contribution around the vicinity of the source hints at a 
stronger dependence of n. In light of these, it is possible that the actual contribution of the interaction term is 
of order nlnn, that is comparable to the term generated as a result of the renormalization process. The detailed 
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analysis of the principal operator, and hence the spectrum, requires developing new approximation methods. 
These questions are postponed to the future works. 
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8 Appendix 



In this section, we will give a brief sketch of the construction of the Lee model and the calculation of the principal 
operator in the light-front coordinate system, and will show that the theory in this coordinate system has the 
same divergence structure. The following oblique coordinate system is chosen, 



u = t + x , 



(139) 



where u is the light-front time coordinate. The infinitesimal invariant distance element, the metric tensor and 
its inverse are also given by 

ds 2 = du 2 - 2dudx - dy 2 - dz 2 , (140) 



9p.v 



I 1 -1 \ 
-10 
0-10 
\ -1/ 



^0-10 Q\ 
-1-10 
0-10 
\ -I) 



The scalar product of the coordinates and the conjugate momenta is 

p M x M = p u u + px + p ± ■ x 1 - 



(141) 



(142) 



where x and x x arc the longitudinal and the transverse coordinates, on the other hand p Ul p and p± are the 
light-front energy, the longitudinal and the transverse momenta, respectively. In the equal-time formulation, 
the bosonic field operator is given by 

The equal-time commutation relations both for fields and for creation and annihilation operators are, respec- 
tively, given by 



[(jKw^x^),^,?/^)] = jsgn(x - yjS^ix 1 - - y^) , 
[o(p, Pi ), aHq, qj.)] = (2n) 3 5(p - q)S^ (p ± - q ± ) . 

The free Hamiltonian of the bosonic sector is 

dp f d 2 p± 



o 2tt J (2tt)2 



■w(p,p_L)a f (p, p±)a(p, p_J . 



(144) 
(145) 

(146) 



where ui{p, p_i_) = — — P 2p v± ■ The positive and the negative frequency parts of the fields evaluated at the point 
zero are given by 



j(+>(0) 
j(->(0) 



dp f d 2 p±_ a(p,p±) 



o 2tt J (2tt)2 ' 
dp f d 2 p ± aJ(p,p±) 



2tt J (2tt) 2 



(147) 
(148) 
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After normal-ordering the creation and annihilation operators, the principal operator takes the form 



®e(E) = f (H - E) u(e) 1 f°° dp f d 2 p ± [°° dq f d 2 q ± 1 1 
A 2 (e) {e) \ A 2 (e) + A 2 (e) 2 J 2tt J (2n) 2 J 2n J (2ir) 2 ^pq a[P ' P±> H - E 



a f (q, qj.) 



Z(e) 



(ffo-B) , Me) 



dp /" d 2 p^ 1 



1 



A 2 (e) A 2 (e) J 2vr J (2tt) 2 2p ff - E + w(p, p ± ) 

dp /" d 2 p ± Z" 00 dg f d 2 q ± 1 + . 1 
S2 I nZ 7^—y^<l±) 



2n J (2tt) 2 J 2tt J (2-k) 2 2^/pq H - E + u(q, q ± ) + u(p, px) 



o(p,pj 



(149) 



We do not need to use any Feynman parametrizations here and only an exponentiation is enough to complete 
the calculations, so the momentum integral in the fourth term in Eq. (|149p is, just, 



dp f d 2 p± 1 



o 2tt J (2tt) 2 2pH - E + lu(p, P± ) 



oo poo 



du 



dp f (Pp± 2u r Ho _ 



o 2^ J (2tt) 



E)p-u(m 2 +p 2 +p 2 L ) 



(150) 



At this stage, we should be careful about the limits of the angular part of the momentum integral. Since we 
work in a coordinate system which covers either the future-cone or the past-cone, after the following change of 
variables, 

p 2 +p 2 ± =s 2 => p = scos9, p ± =ssin6>, (151) 
the integration interval of the ^-integral becomes [0, -|]. Equation (|150|) is, then, 



1 



(2tt) 3 



rr/2 

dss 2 I d9 sin 9 / d(f>e 



8(2^) 2 



JO 

ra 2 u poo 



-2u(H -E)s cos 6~us 2 



du ■ 



1 



dsse- s Z 4 -^ 
o ^I(H - E) 



1 - e - s V^(Ho-E) 



(152) 



By using the exponential representation of the fractions in the fifth term in Eq. (|149p , the principal operator is 
given by 



A 2 (e) 



Z(e). 



(H -E) , Me) 



1 



A 2 (e) A 2 (e) 32tt 

oo />oo 



du — jj^- I dsse~ sJ/4 



1 _ e -ST/S(H -E) 



V^(Ho - E) 



• I • r a* r dp f d2p± r dq f d2q± 



x e 



-qa 2 -soj(q,q ± ) -p/3 2 -su>(p,p ± ) t 



a T (g,qx)e 



-s(H -E) 



a(p,P±) 



(153) 



With the help of the redefinitions of the mass and the coupling constant below 



Mf) = A*fl 1 
A 2 (e) A| 32^ 2 

1 1 1 



-urn r oo 



dlt ■ 



ds S 2 e- s2 / 4 . 



A 2 (e) A| 64^ 2 



-ttm i-oo 



du • 



(154) 
(155) 



one can take the limit e — ► + after dividing both sides by Z(e) and hence the renormalized principal operator 
takes the form, 



ggCg) _ (Hp - £?) u fl _J_ 

A 2 fl 32^ 2 y 



4 



-m u r°° 



du ■ 



dsse- s /4 - 

V^(H -E) 



1 - sVM# - -B) + ^V#o - E) — e -^{H -E) 



2 

x e 



OO /"OO 



d./ da / d(3 I ± f ^ r dq f ^ 



2tt J (2tt) 2 J 2tt 7 (2?r) 



-<ja 2 -sw(g,q i ) -p/3 2 -sw(p,p ± ) f 



a T (<?, qj_)e" 



-s(Ho-E) 



a(p,P±) ■ 



(156) 
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Now to see the divergence patterns, we can again calculate the bare mass and the bare coupling constant 
asymptotically in e, as a result we find the following, 

- " 0+ ' <157) 

W)' J k + ^' ne as E " 0+ ' (158) 

We note that the divergences are controlled by the cut-off parameters in exactly the same way as in the previous 
cases. We believe this is in a cord with the discussion presented by the authors in 25j about the equivalence 
of the covariant perturbation theory and the light-front perturbation theory. This may be seen as another 
verification of this equivalence at a nonperturbative level. 

The asymptotic limit of the renormalized principal operator can, of course, be analyzed in this case, as well. 
For the calculations to be done for this analysis repeat themselves, we will not continue further in this direction. 
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